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SUMMARY 

A successful design of a nuclear aircraft propulsion system must 
meet all specifications of performance and safety at a minimum weight. 
Heat pipes utilized for removing core heat may be able to reduce sys- 
tem weight by eliminating the need for heavy pumps and high pressure 
piping. 

A contract was awarded to C. Silverstein for the study of heat pipe 
applications to nuclear aircraft propulsion systems (ref. 1). This study 
was limited to the thermal and mechanical design of heat pipes and did 
not include core criticality calculations. 

Therefore, a study was conducted by Lewis Research Center con- 
currently to perform the neutronic calculation on the proposed heat pipe 
reactor design, The study conducted revealed that utilizing heat pipes for 
a nuclear airplane reactor application appeared promising when heat 
pipe performance was applied to the limit of heat pipe technology. The 
design parameters calculated are 
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Number of heat pipes , 2740/half core 

Radial heat flux 0. 88xl0 6 Btu/hr-ft 2 

(2. 77x10® watts/meter 2 ) 

Core diameter 45 inches (114 cm) 

[0.075 in. (0. 18 cm) fuel,, 0. 035 in. (0. 089 cm) clad] 

Core L/D . 1.0 

Heat pipe vapor temperature 2100° F (1422 K) 

Fuel enrichment 37. 5 percent 

Fuel loading 3770 kg 

Clad temperature 2300° F (1530 K) 

Clad stress 10 000 psi (6890 N/cm 2 ) 

Core K e ff . 1.1 

Power 300 MW 


The active core diameter of 45 inches (114 cm) can be decreased 
by increasing the enrichment and reducing the fuel loading of the core. 
This would be desirable since control drum application appears difficult 
at 45 inch (114 cm) outside diameter due to the reduced boundary flux. 


INTRODUCTION 


Many studies have been conducted where heat pipes are applied to 
nuclear reactor cores or thermionic energy conversion devices. These 
studies indicate possible advantages in efficiency and design simplification. 
Weight reduction has also been accomplished through the elimination of 
pumps and high pressure piping. 

On January 22, 1969 a contract was awarded to Mr. C. C. Silverstein 
for the study of heat pipe applications in nuclear aircraft propulsion 
systems. This study included a preliminary design evaluation of a heat 
transport system in which heat pipes were used for removing and trans - 
porting the reactor core heat to the aircraft engine air heaters. This 
study was limited to the thermal and mechanical design of heat pipes and 
did not include core criticality calculations. 
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Therefore, a reactor design study was conducted by Lewis Research 
Center concurrently with Mr. Silverstein. Its purpose was to provide 
reactor neutronic information which could be used, together with heat 
pipe information, to establish reactor size and weight and subsequent 
shielding size and weight. The neutronic calculations determine the core 
size and fuel enrichment for a core of 1,1, given the heat pipe 

design characteristics. This report presents the results of this study. 

SYMBOLS 

Aq internal heat generation rate, Btu/hr-cu. in. (watt/cc.) 

Cj, C£ integration constant 

D core diameter, in. (cm) 

eV neutron energy, electron volts 

K thermal conductivity, Btu/hr-ft-°F (watts/meter-K) 

n number of heat pipes 

L heat pipe length, core length, in. (cm) 

q heat flux, Btu/hr-ft^ (watts/meter^) 

r radius, in. (cm) 

S r radial spacing, in. (cm) 

t thickness, in. (cm) 

t temperature, °F (K) 

u neutron energy, lethargy 

Subscripts: 

A Axial 

cell heat pipe fuel cell 

cl clad 

core active reactor core 
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f fuel 

r radial 

v vapor 

GENERAL DESCRIPTION OF COMPONENTS 

A sketch of the proposed design is presented in figure 1. The core 
is cylindrical. It consists of heat pipe fuel elements surrounded by 
four inch (10. 2 cm) thick side and end reflectors. The heat pipes start 
at the core midplane, exit through the top and bottom end reflectors and 
terminate in heat exchangers just outside the end reflectors. 

The heat pipes consist of a fluid, wick and tube. The part of the heat 
pipe within the core is surrounded by an annulus of fuel which in turn is 
surrounded by cladding, This fueled portion covers about half of the 
active length of the heat pipe. The heat pipe fuel elements are clustered 
tightly such that nominally zero spacing occurs between their outer 
diameters within the active core region. With the fuel and cladding 
terminating at the end reflectors, there is space between the heat pipes 
in the heat exchangers. The heat pipes are supported by the end reflectors. 
These reflectors are in turn supported by a containment vessel. 

The heat pipes and reactor core were designed so that a minimum 
amount of design feedback was involved. For example, the heat pipes 
were surrounded by the fuel rather than vice versa so that once the heat 
pipes were designed a change in fuel thickness did not necessitate recal- 
culation of the heat pipe. 

HEAT PIPE AND REACTOR DESIGN CHARACTERISTICS 

The initial heat pipe design proposed by Silverstein resulted in a 
much too large active core diameter. A second heat pipe design was 
satisfactory. These heat pipes are circular in cross section and consist 
of a liquid metal vapor space, a capillary wick, and a liquid annulus all 
enclosed in a tube (fig„ 1). They are required to remove 300 MW of 
reactor thermal power at a maximum vapor temperature of 2100° F (1422 K). 
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The heat pipe design is reported in reference 1. A summary of the 
design is as follows: 


Materials . . . . , 

Heat pipe cross section . . . 


Heat pipe length. . 

Average temperatures .... 

Maximum axial heat flux . . . 

Average axial heat flux . .. .. 

Average heat transport . . . 
rate/heat pipe 


Heat pipe fluid - Sodium 
Wick and wall - TZM moly 

Circular consisting of inner vapor 
space of radius 0. 2 5 inches (0, 636 
cm) surrounded by annular TZM 
wick 5 mills (0. 013 cm) thick with 
30 percent void and effective pore 
diameter of 2 microns. Wfck is 
surrounded by 10 mils (0.025 cm) 
thick annulus filled with sodium 
(which also occupies wick void 
volume). The sodium annulus is 
surrounded by TZM containment 
wall with thickness of 33 mils 
(0. 084 cm). 

20 Inch (50. 8 cm) evaporator length 
in each half of core 

Heat pipe vapor - 2000° F (1367 K) 

307 kW/in. 2 (47. 6 kW/cm 2 ) 

279 kW/in. 2 (43. 3 kW/cm 2 ) 

54. 9 kW 


Surrounding each heat pipe is a fuel annulus, clad on its outside 
diameter. Both the thickness of the fuel and its enrichment can be 
varied. For calculating the required fuel thickness, fuel enrichment, 
and core size the following reactor design criteria was established: 
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Core fuel UN 

Clad TZM 

Core K eff . . . 1.1 

Core lifetime 10 000 hr 

Peak fuel burnup 5 percent 

Fuel thickness 0. 075 in. (0. 19 cm) 

Clad thickness . 0. 035 in. (0. 089 cm) 

Active cpre diameter . 45 in. (114 cm) max 

Side and end reflector 4 in. (10. 2 cm) - Moly 

Radial peak to average power ratio 1.1 

Axial peak to average power ratio 1.0 


DESIGN PROCEDURE 

Using the heat pipe design outlined above and the reactor specifi- 
cations, the number of heat pipes, radial heat flux, core diameter, 
heat pipe fuel cell radial temperature gradient, fuel loading, and 
enrichment can be calculated. The following is a presentation of the 
calculations of these parameters. 

Number of Heat Pipes 

With the heat pipes removing core heat through both end reflectors, 
the reactor can be designed considering symmetry about its midplane. 
Therefore, 150 MW of heat must be removed from each half of the core. 

The required number of heat pipes for one half core is expressed 
as follows: 

150 MW x 4 
«A* D v 


n = 


(i) 
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where D is the vapor diameter of the heat pipe. For a vapor diameter 

v 2 
of 0. 5 inch (1. 27 cm) and an average axial heat flux of 279 kW/in. , then 

2740 heat pipes are required in each half the core. 

Core Diameter 

The active core diameter is dependent upon the number of heat pipes, 
fuel thickness and cladding thickness. Since the fuel and cladding thick- 
ness have been specified the core diameter has in effect been dictated by 
the 2740 required tubes. 

Figure 4 presents the proposed design of the heat pipe fuel cell. The 
area of this hexagonal cell may be calculated as follows: 



Multiplying equation (2) by the number of heat pipes gives the active core 
area if the scallops of the outer perimeter are neglected, 

A core sn Y( 2r 6 + S r) 2 <» 

From equation (3) the diameter of the core is 

D 

core 


( 4 ) 

or in terms of fuel thickness tj and cladding thickness t j, equation (4) 
becomes 
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Finally, for an assumed of 0.035 inch (0.089 cm), an assumed 
of 0.075 inch (0. 19 cm), a S r of 0, and 2740 tubes, the core diameter 
is 45 inches (114 cm). This value meets the maximum core diameter 
design specification^ A smaller core diameter could be obtained by re- 
ducing the clad and/or fuel thickness or by making the fuel cross sec- 
tional geometry hexagonal instead of circular. 

Radial Heat Flux 

The heat removed axially from the core enters the heat pipes ra- 
dially through its walls. The radial heat flux is dependent upon the axial 
heat flux and the pipe’s length and diameter or 

CLA^l ^A r l 

q = JL_! = JL! (6) 

r 2 tt^L 2L 

The reactor core requirement of a L/D of 1 has set the heat pipe 
length for half a core at approximately 20 inches (50. 8 cm). For a q^ 
of 279 kW/in.^ (43. 3 kW/cm^) of vapor cross section, the q r is 0. 88*10® 
Btu/hr-ft^ (2.77x10® W/m^). The radial heat flux directly effects both 
the maximum fuel temperature and thermal stresses of the fuel. The se- 
lection of UN as the fuel reduces this effect over fuels such as UOg due 
to UN’s higher thermal conductivity. 

Heat Pipe Fuel Cell Radial Temperature Gradient 

The axial heat removed by the heat pipes is generated in the UN fuel 
surrounding the heat pipe tube wall. This internal heat generation and 
resultant radial flux results in a radial temperature gradient from the 
heat pipe liquid and wall interface through the tube wall, fuel, and to the 
clad wall interface around the fuel. No temperature gradient exists in 
the cladding since no heat is transferred. 
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The total temperature gradient from the fuel outer wall to the heat 
pipe inner wall is 



(See Appendix A for derivation. ) 

A radial heat flux of 0. 88 Btu/hr-ft 2 (2. 77x10^ W/m 2 ) results in a 
AT of 200°. For a maximum heat pipe vapor temperature of 2100° F 
(1422 K) the resultant clad temperature is 2300° F (1530 K). 

For TZM at 2300° F (1530 K), the 10 000 hour creep-rupture 
stress is about 10 000 psi (6890 N/cm^) f Since the clad has no thermal 
gradients (no heat is lost radially outward) its stresses would be limited 
to that caused by fuel swelling, fuel growth, fission product gases, etc. 

Fuel Loading 

The 0.075 inch (0. 19 cm) thick fuel surrounding each heat pipe also 
sets the total fuel loading of the core in kg of weight. With the reactor 
core required to operate for 10 000 hours at 300 }MW, the amount of fuel 
that will be used can be calculated as follows: 

300X10 6 Watts x 10 000 hr x 3600 sec/hr = 10. 8xl0 15 Watt-sec 

For U 235 fuel the burn-up rate is 1.452x10 ** gm U 23 ^/Watt-sec of op- 
eration. The total amount of fuel used is 

10.8X10 15 Watt-sec x 1.452X10 -11 gm U 235 /Watt-sec = 157 kg 

From figure 7 a 0. 075 inch (0. 19 cm) thick fuel for 5480 tubes (total 
core) results in a core loading of 3770 kg. This loading would result in 
a peak burn-up of 4. 2 percent which is less than the design goal of 5 per- 
cent. 
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Now all thermal and physical design variables have met the require- 
ments of the section titled ’’Heat Pipe and Reactor Design Character- 
istics. ” The final design variable to satisfy is the enrichment for a core 
K e ff = !• 


CORE CRITICALITY CALCULATIONS 

The code used in the neutronic analysis of this report is a two- 
dimensional discrete angular segmentation transport program referred 
to as TDSN (ref. 2) . It is a numerical iterative finite difference method 
in which the continuous angular distribution of neutron velocities is rep- 
resented by considering discrete angular directions. The output of the 
transport program includes the core multiplication factor, integrated 
power ratio and flux shapes for each energy group. 

A required input to the TDSN transport program are the core cross 
sections. Multigroup cross sections used in this analysis were obtained 
from GAM II and GATHER- II programs (refs 3 and 4). These are then 
converted to macroscopic cross sections by a Lewis Research Center 
written code entitled MACROS. 

For the calculations of this report the core was divided into one ra- 
dial zone and one axial zone. Seven group microscopic cross sections 
were obtained from the GAM II and GATHER II programs over the energy 
ranges presented in table I. The macroscopic cross sections calculated 
by MACROS used the atom densities listed in figure 4. A Pq approxi- 
mation for neutron scattering and a Sg discrete angle approximation 
were used in the neutronic calculations. The core K e ^ was calculated 
over a large range of enrichments. For each fuel enrichment the core 
atom densities were readjusted. The results of these calculations are 
presented in figure 3, core as a function of Percent Fuel Enrich- 

ment. The required core K e ^ was achieved at a 37. 5 percent fuel en- 
richment. 

The normalized power distribution is plotted in figure 4. The plot 
shows that the radial peak to average power ratio is approximately 1.45 
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compared to a desired 1,1. Since the fuel enrichment is low, however, 
fuel zoning could be applied to better meet this requirement. The low 
flux at the boundary also results in a low radial neutron leakage of 10 per- 
cent. This would negate the use of boundary control such as control 
drums. However, if the radial power distribution were flattened to 
nearer the 1. 1 ratio the radial leakage would increase and boundary con- 
trol may then be feasible. 

SUMMARY OF RESULTS 

The values of the design variables that meet the requirements set 
forth in the section titled "Heat Pipe and Reactor Design Characteristics” 
section are 

Number of heat pipes . 2740/half core 

Radial heat flux 0. 88xl0 6 Btu/hr-ft 2 

(2.77x10^ watts/meter 2 ) 

Core diameter .............. 45 inches (114 cm) 

(0.075 in. (0. 19 cm) fuel, 0, 035 in. (0.089 cm) clad^ 

Core L/D ....... ....... .......... 1.0 

Heat pipe vapor temperature . 2100° F (1422 K) 

Fuel enrichment . 37. 5 percent 

Fuel loading 3770 kg 

Clad temperature 2300° F (1530 K) 

Clad stress ............ 10 000 psi (6890 N/cm 2 ) 

Core o ............. . 1.1 

Power ........................ 300 MW 

The heat pipe design may be compared with a liquid metal fast reac- 
tor design (ref. 5) and a water moderated thermal reactor design (ref, 6). 
The comparison is qualitative. The reactors were designed for the same 
power and life but for different burnups. The 45-inch (114-cm) diameter 
at a 4. 2 percent burnup for the heat pipe core compares with a 30- inch 
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(76.2-cm) diameter fast reactor at a 12 percent burnup and 62-inch 
(157-cm) diameter thermal reactor at a 20 percent burnup. If the burn- 
up rate for the heat pipe design were increased to 10 percent, the core 
diameter reduces to 39 inches (99 cm). This may be possible since a 
fuel enrichment of 37, 5 percent could be increased to maintain a core 

K eff of lo 

The scope of this study did not include the design of control rods or 
drums. The calculated boundary neutron flux and radial leakage, how- 
ever, indicate that control drums would be difficult to apply. Therefore, 
a more optimized design should endeavor to flatten the radial power dis- 
tribution, reduce the core diameter, or incorporate poison heat pipe 
tubes within the core. 


CONCLUSIONS 

The following conclusions are made based upon the preliminary 
analysis conducted in this report. 

1. The concept of utilizing heat pipes for the removal of core heat in 
a nuclear aircraft 300 MW reactor appears promising when the heat pipe 
performance is pushed to the limit of present heat pipe technology. 

2. The active core diameter of 45 inches (114 cm) (4, 2 percent burn- 
up) compares with 30 inches (76. 2 cm) for a liquid metal fast reactor 

(12 percent burnup) and 62 inches (157 cm) for a water moderated thermal 
reactor (20 percent burnup). A smaller core diameter of 39 inches (99 cm) 
is possible for the heat pipe core if a 10 percent burnup were assumed. 

3. Control drum application for the current design will be difficult due 
to the reduced boundary flux and high energy level of the neutrons. 

4. Further reduction in fuel loading and subsequent core diameter is 
possible by reducing the 0.075 inch (0. 19 cm) fuel thickness and the 0.035 
inch (0. 089 cm) clad thickness. 
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APPENDIX A 


For the case of heat production in a cylinder 


‘IriUa 


r dr V dr/ K 


dfr*V-— ’ 

\ dr / K 


ctt = _ + ^1 

dr 2K r 


The boundary conditions for equation (Al) are 


— - 0 at r = r. 


c i = 


A r^ 
A 0 r 5 


dt_ _ _;^V [ 1 A Q r 5 
dr 2K r 2K 


integrating equation (A2) gives 


t = —/-— + rj? In r\ + C ? 
2K 1 2 5 ( 2 


The boundary conditions for equation (A3) are 
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Equation (A4) represents the temperature drop from the heat pipe tube- 
outer wall radially outward to the fuel outer wall. For the total tempera- 
ture gradient from the fuel outer wall to the heat pipe inner wall (ref. r^ 
and r 3 of fig. 2). 



Substituting equation (A5) into equation (A4) and applying the boundary 
cpndition t = tg at r 5= r g gives 
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TABLE I, - NEUTRON ENERGY GROUPS 


Group 

Neutron energy range 


eV 

u 

1 

1.49X10 7 to 2, 23xl0 6 

-4, 00X10" 1 to 1.5 

2 

2. 23X10 6 to 8.2X10 5 

1. 5 to 2. 5 

3 

8.21X10 5 to 7.10X10 3 

2.5 to 7.25 

4 

7, 10X10 5 to 7.49X10 2 

7,25 to 9. 50 

5 

7.49X10 2 to 29, 0 

9. 50 to 12. 80 

6 

29.0 to 0.414 

12,80 to 17.0 

7 

0.414 to 0 

17.0 to — 
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